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Q-l Abstract 

<D . 

r-j _ We construct a new example of the spinning-particle model without 

Grassmann variables. The spin degrees of freedom are described on the 
base of an inner anti-de Sitter space. This produces both T M and r MI/ - 
matrices in the course of quantization. Canonical quantization of the 
model implies the Dirac equation. We present the detailed analysis of 
both the Lagrangian and the Hamiltonian formulations of the model and 
obtain the general solution to the classical equations of motion. Compar- 
ing Zitterbewegung of the spatial coordinate with the evolution of spin, we 
ask on the possibility of space-time interpretation for the inner spin-space. 
We enumerate similarities between our analogous model of the Dirac equa- 
tion and the two-body system subject to confining potential which admits 
only the elliptic orbits of the order of de Broglie wave-length. The Dirac 
equation dictates the perpendicularity of the elliptic orbits to the direction 
of center-of-mass motion. 
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1 Introduction 

Although the true understanding of spin is achieved in the frame- 
work of quantum electrodynamics, a lot of efforts has been spent in 
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attempts to construct the relativistic mechanical model of a spin- 
ning electron [2-19]. Just after the introduction of the electron spin 
as a new quantum number by Pauli pQ, Uhlenbeck and Goudsmit 
suggested its naive interpretation in terms of the inner angular mo- 
mentum P]. The first papers devoted the semiclassical description 
of spin dynamics can be traced back up to FrenkeQ |4j and Thomas 
[5]. Bargmann, Michel and Telegdi have demonstrated [7] that the 
models practically exactly reproduces the spin dynamics of polarized 
beams in uniform fields. However, the miracle is that the models 
based on these schemes do not produce the Dirac equation through 
the canonical quantization. One possible solution to the problem has 
been found by Berezin and Marinov [H [9] in the framework of Grass- 
mann mechanics. The problem here is that the Grassmann mechan- 
ics represents a rather formal mathematical construction. It leads 
to certain difficulties [HI EI] in attempts to use it for the description 
of spin effects on the semiclassical level, before the quantization. We 
also point out that there is no generalization of Grassmann mechan- 
ics on higher spins [TO] . 

Hence it would be interesting to construct the non-Grassmann 
model of the Dirac equation. While the problem has a long history 
(see [2-7, 10-19, 23, 24, 26] and references therein), there appears 
to be no wholly satisfactory solution to date. Our believe on the 
existence of such a kind model is based on the following well-known 
observation. The Dirac spinor \1/ can be used to construct the four- 
dimensional current vector, ^r^ 1 !/, which is preserved for solutions 
to the Dirac equation, d^^T^^) = 0. Hence its null- component, 
ty'ty > 0, admits the probabilistic interpretation, and we expect 
that one-particle sector of the Dirac equation admits description in 
the framework of relativistic quantum mechanics (RQM). So we can 
look for the corresponding semiclassical model, which would lead to 
the Dirac equation in the course of canonical quantization. 

However, it is well-known that adopting the RQM interpretation, 
we arrive at the rather strange and controversial picture. To remind 
this, we use the Dirac matrices a % and /3, to represent the Dirac 
equation in the form of the Schrodinger equation 

ihd t ^ = H^, H = ca% + mc 2 (3. (1) 
Then H may be interpreted as the Hamiltonian. If we pass from the 

1 See a modern revision of Frenkel's paper in [6]. 
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Schrodinger to Heisenberg picture, time derivative of an operator is 
iha = [a, H}. For the basic operators of the Dirac theory we obtain 

±i = can, thai = 2(cpi - Hcti), pi = 0. (2) 

Some immediate consequences of these equations are enumerated 
below. 

• The wrong balance of the number of degrees of freedom. Ac- 
cording to the first equation from (jSJ), the operator ca % rep- 
resents velocity of the "center of charge" x % [11) . Then the 
physical meaning of the operator p % became rather obscure in 
both the semiclassical and the RQM framework. 

Various approaches to the problem has been considered in the 
literature. Schrodinger noticed [UJ that besides the center of 
charge, x, in the Dirac theory we can construct the "center- 
of-mass" operator = X{ + \%hcH~ x OLi in such a way, that 
p % turns out to be the mechanical momentum for x, x ~ p 
(various versions of this operator have been discussed in the 
works [TBI [Til [T5]). Following this way, Schrodinger assumed 
the naive interpretation of the Dirac electron as a kind of com- 
posed system (we return to this subject in Subsection 6.1). In 
contrast, Foldy and Wouthuysen [H] assumed that x l does not 
correspond to an observable quantity. 

• Free electron follows complicated trembling trajectory. The 
equations (j2J) can be solved, with the result for x l (t) being 

[rung 

x* = cj + b P H + c i exp(-'^-t). (3) 

The trajectory is a superposition of the rectilinear motion along 
a straight line, Oj + bprf, and the rapid oscillations with higher 
frequency W- ~ The oscillator motion is called Zitterbe- 

wegung. Schrodinger had compared spin with the angular mo- 
mentum associated to the Zitterbewegung. He found [UJ that 
they differ by the factor 2 and concluded that spin can not 
be identified with the Zitterbewegung. In our model, we would 
be able to construct the variables for which the identification 
turns out to be possible, see Subsection 6.1. 
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• Since the velocity operator ca l has eigenvalues ±c, we con- 
clude that a measurement of a component of the velocity of a 
free electron is certain to lead to the result ±c. Besides, since 
the operators a 1 do not commute, components of velocity in 
different directions can not be instantaneously measured. 

In view of this, one may assume that the Dirac equation does not 
admit the RQM interpretation. Other possibility might be that the 
basic operator x appeared in the Dirac equation do not correspond 
to the physically observable quantity. This possibility is supported 
by the seminal work of Foldy and Wouthuysen [H] , where they have 
constructed (in a Lorentz non-covariant manner) the position oper- 
ator with reasonable properties. In particular, it obeys the equation 
X = This leads to further complications, as the Dirac equation 
gives no evidence which of these two operators should be identified 
with the position of an electron. We return to the subject in Sec- 
tion 5, where we propose the Lorentz-covariant classical- mechanical 
analog for the Foldy- Wouthuysen operator. 

To understand the controversial properties of the one-particle 
Dirac equation, it would be desirable to have at our disposal the 
spinning-particle model which leads to the Dirac equation in the 
course of canonical quantization. We construct and discuss the 
analogous model of such a kind in this work. It shows the same 
properties as those of the Dirac equation in the RQM interpreta- 
tion. Analyzing the present model, we have been able to identify 
the origin of the problems, see Section 7. The modified model which 
turns out to be free of the problems mentioned above has been pro- 
posed in the recent work [23] . 

The work is organized as follows. The operators of the Dirac 
theory which are associated with spin are the r M -matrices as well 
as the Lorentz generators T^ u . Their commutators can be identified 
with the Poisson-brackets of angular momentum of five- dimensional 
space with the metric r\ = (— ,+,+,+,— ) [23]. So the space can 
be taken as the underlying configuration space for the description 
of spin. As the number of independent angular-momentum compo- 
nents is less than dimension of the spin phase space, dynamics of 
spin should be restricted to an appropriate subspace which we call 
the spin surfaced f23]. The corresponding configuration space turns 

2 For the case of non-relativistic spin | 24 |, the surface has natural structure of the Hopf 
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out to be the anti-de Sitter space. To make the work self-consistent, 
we present the construction of the spin surface in Section 2. 

Hamiltonian formulation of the model with spin sector of such 
a kind has been announced in [25]. In the Sections 3, 4 and 5 
we carry out the detailed analysis of Hamiltonian and Lagrangian 
formulations of the model and analyze it on the classical as well as 
on the quantum level. In the Section 6 the classical equations of 
motion are integrated and discussed in details. Section 7 is left for 
the conclusions. 

2 Algebraic construction of relativistic spin sur- 
face 

We start from the Dirac equation written in the manifestly- 
covariant form 

(p^ + mc)*^) = 0, (4) 

where p^ = —ihd^. We use the representation with hermitian T 
and antihermitian P 

then [T^,T U ] + = -Tafi = (~ + ++)> and r ° r N r ° are the 
Dirac matrices a\ (3 [TJ]. We take the classical counterparts of the 
operators and p^ = —iTid^ in the standard way, which are x M , p u , 
with the Poisson brackets {x fl ,p u } = r]^. 

Let us look for the classical variables that could produce the 
T -matrices. According to the canonical quantization paradigm, the 
classical variables, say z a , corresponding to the Hermitian operators 
z a should obey the quantization rule 

[z a ,zP] = ih{z<*,z%^. (6) 

In this equation, [ , ] is the commutator of the operators and { , } 
stands for the classical bracket!!. To avoid the operator-ordering 
problems, we will consider only the sets of operators which form 

fibration S 3 -> S 2 , see |25j . We are grateful to Prof. A. Nersessian for the comments on this 
subject. 

3 It is the Poisson (Dirac) bracket in a theory without (with) second-class constraints. 
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the Lie algebra, [z a , z 13 } = c a/3 7 5 7 . So our first task is to study the 
algebra of T -matrices. We note that commutators of T M do not form 
closed Lie algebra, but produce 50(1,3) Lorentz generators 

[T", V] = -2iF^, (7) 

where = ifP-PT"). The set V 1 , forms a closed algebra. 
Besides the commutator (J7J), we have 

[r^,T a ] = 2i{rf a T v -r)™^), 
[P*", T al3 } = 2i(if c T vf> - r]^T va - rf a V^ + ^P**). (8) 

The algebra can be identified with 5*0(2,3) Lorentz algebra with 
generators J AB : 

[J AB , J CD ] = 2i( V Ac j BD - rj AD J BC - 

v BCjAD + J] BDjAC^ (g) 

assuming = J 5 ^, P" = 

To reach the algebra starting from a classical-mechanics model, 
we introduce ten-dimensional "phase" space of the spin degrees of 
freedom, u A , n B , equipped with the Poisson bracket {u a ,tc b } = 
r] AB . Consider the inner angular momentum 

J AB = 2(U A 7T B -U B 7T A ). (10) 

Poisson brackets of these quantities form the algebra 

{J AB , j cd } pb = 2( v ac j bd - v ad j bc - 

r] BCjAD + 7] BD J AC ) _ (n) 

Comparing (fTTj) with (jUJ) we conclude that the operators T M , 
could be obtained by quantization of J AB . 

Since J AB are the variables which we are interested in, we try 
to take them as coordinates of the space u A ,n B . The Jacobian 
of the transformation (u A , tt b ) — > J AB has rank equal seve rfl. So, 
only seven among ten functions J AB (u, ir), A < B, are independent 
quantities. They can be separated as follows. By construction, the 
quantities ffTUl) obey the identity e^ val3 J 5 U J a p = 0, this can be solved 
as 

jij = (j50yl(j5ij0j _ J5jj0iy ( 12 ) 



The rank has been computed using the program: Wolfram Mathematica 8. 
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Hence we can take J 5m , J° l as the independent variables. We could 
complete the set up to a base of the phase space (u A , n B ) adding 
three more coordinates, for instance u 3 , u 5 , 7r 5 . Quantizing the 
complete set we obtain, besides the desired operators J 5m , J°\ some 
extra operators a) 3 , a) 5 , 7r 5 . They are not present in the Dirac the- 
ory, and are not necessary for description of spin. So we need to 
reduce the dimension of our space from ten to seven imposing three 
constraints which we denote T a (o;,7r) = 0, a = 3,4,5. We require 
the surface defined by the constraints be invariant under action of 
the angular-momentum generators, that is 

{T a (co,7r),J AB } = 0. (13) 

The only quadratic SO (2, 3) -invariants which can be constructed 
from co A , ii B are u a u>a, uj a ita and t^ a t^a- So we restrict our model 
to live on the surface defined by the equations 

T 3 = n A 7i A + a 3 = 0; (14) 



T 4 ee u a ua + a 4 = 0, T 5 ee uj a tta = 0, (15) 

where 0,3, a 4 are positive number^]. It is called the spin surface. 

The first equation from ffl5|) states that the configuration space 
of the spin degrees of freedom is anti-de Sitter space. 

In the Hamiltonian formulation, the equations T a = appear 
as the Dirac constraints. So, we classify them in accordance with 
their algebraic properties with respect to the Poisson bracket. The 
brackets read 

{T 3 , T 4 } = -4T 5 , {T 3 , T 5 } = -2T 3 + 2a 3 , 

{T 4 ,T 5 } = 2T 4 -2a 4 . (16) 



Taking the combination 



f 3 ee T 3 + ^T 4 , (17) 
a 4 



we have the algebra 



{f 3 , T 4 } = -4T 5 , {f 3 , T 5 } = -2T 3 + 2^T 4 , 

a 4 



5 Tho positivity guarantees the Causal dynamics of our particle, see Eq. I I72II . 
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{T 4 ,T 5 } = -2a 4 + 2T 4 , (18) 

that is Poisson bracket of T 3 with any constraint vanishes on the 
surface, while the Poisson brackets of constraints T 4 , T 5 form an 
invertible matrix on the surface. In the Dirac terminology, the set 
T 3 is the first-class constraint, while T 4 , T 5 is a pair of second-class 
constraints. 

There are various reasons to take the functions T a invariant. 

A) Consistency of canonical quantization of a system with second- 
class constraints implies replacement the Poisson by the Dirac bracket, 
the latter is constructed with help of the constraints. For the mo- 
mentum generators it reads (i = 4, 5) 

{J AB , J CD } DB = {J AB ,J CD } - {J^HT,,!}}- 1 ^, J CD }. (19) 

If the surface is invariant, the second term on the r. h. s. vanishes, 
and the Dirac bracket of J AB coincides with the Poisson bracket, 
Eq. (ITTj) . So, as before, we have the desired algebra. 

B) Presence the first-class constraint T 3 implies that we deal with a 
theory with local symmetry. Generators of the symmetry are pro- 
portional to the constraints [211 123 E2] • Suppose that the first-class 
constraint is not invariant, that is we have {T 3 , J ab }pb 7^ 0. It 
should imply that the variables J AB are not inert under the local 
symmetry, 5J AB ~ {T3, J ab }pb 7^ 0. Hence the non-invariant con- 
straints would be responsible for gauging out some of the variables 
J AB , which is not under our interest now. 

Let us discuss convenient parametrization of the spin surface. 
The matrix — ^q^a^b^'^ - has rank equal ten. So the quantities 

J 5 », J°\ T 4 , T 5 , u 5 , (20) 

can be taken as coordinates of the space (u A , n B ). The equation 
jab _ 2( y co A n B - cu B n A ) implies the identity 

J AB J AB = S[{u A f{n B f - (A) 2 ] = 

8[(T 4 -a 4 )(T 3 -a 3 )-(T 5 ) 2 ], (21) 



then the constraint T 3 can be written in the coordinates (|20|) as 
follows: 

8(i 4 — a 4 J 
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where J Jjf are given by Eq. ( !T2|) . Note that T 3 does not depend on 
u 5 . On the hyperplane T 4 = T 5 = it reduces to 

- 8a 4 T 3 = ( J AB f - 8a 3 a 4 = 0. (23) 

Eq. (12"3"1) states that the value of £0(2, 3) -Casimir operator (J AB ) 2 
is equal to 80304. In quantum theory, for the operators we 
have: J AB Jab — 20o 2 . So we can fix the product of our parameters 
as 

a 3 «4 = — • (24) 

As we have pointed out above, the function T 3 represents gener- 
ator of local symmetry. The coordinate u 5 is not inert under the 
symmetry, 5u 5 ~ {T 3 ,c<; 5 } 7^ 0. Hence cu 5 is gauge non-invariant 
(hence non-observable) variable. 

Summing up, we have restricted dynamics of spin on the sur- 
face (P3J, (TT5|) . If (12"U|) are taken as coordinates of the phase space, 
the surface is the hyperplane T 4 = T 5 = with the coordinates 
J 5fl ,J \u 5 subject to the condition ff23l . Since a; 5 is gauge non- 
invariant coordinate, we can discard it. It implies that we can quan- 
tize J 5m , J°* instead of the initial variables u A , n B . 

Following the canonical quantization paradigm, the variables must 
be replaced by Hermitian operators^] with commutators resembling 
the Poisson bracket 

U=tt{, }lj-J- (25) 

Similarly to the case of T -matrices, brackets of the variables J 5m , 
J 0t do not form closed Lie algebra. The non closed brackets are 

{ J 5i , J^} = {J°\ J°i} = -2jv, (26) 

where J u is given by Eq. ( TT2|) . Adding them to the initial variables, 
we obtain the set J = (J 5m , J°\ J u ) which obeys the desired 
algebra (ITTj) . 

According to Eqs. (El), (TTTT) the quantization is achieved replacing 
the classical variables J 5fl , J^ u on r-matricefl We assume that 

6 The matrices T^, r M " are Hermitian operators with respect to the scalar product 

(*i,* 2 ) = *Ir°* 2 . 

7 Replacing i l 1 2[ ) by an operator J IJ (IV , F 0t ) we arrange the operators T in such a way, that 
jij (p) _ pij i 
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u A has a dimension of length, then J AB has the dimension of the 
Planck's constant. Hence the quantization rule is 

This implies that the Dirac equation (TJJ can be produced by the 
constraint 

T 2 = V(l + mch = 0. (28) 



3 Hamiltonian formulation 

Our next task is to formulate the variational problem for our 
spinning particle. As it has been discussed above, we need a theory 
which implies the constraints (JHJ), ({15]) and (|28|) . Since they are 
written on the phase space, it is natural to start from construction 
of an action functional in the Hamiltonian formalism. 

Hamiltonian action of non-singular (that is non-constrained) sys- 
tem reads J dr[PQ — H(Q, P)]. When the theory is singular, Hamil- 
tonian action acquires more complicated form [2U1 ETj 122]- To re- 
mind its structure, as well as to justify our choice of the action 
fl38|) . f )39|) . we outline the Hamiltonian formulation for singular La- 
grangian theory of a special form. For our purposes it will be suffi- 
cient to discuss the Lagrangian 

S = J drL(Q, Q, e a ), (29) 

that is the configuration variables are divided on two groups, Q 
and e a , where e a enter into the action without derivatives. We also 
suppose that 

d 2 L 

det — , — r- ^ 0. (30) 
dQdQ r 

Following the standard prescription, we construct Hamiltonian for- 
mulation for the action (|2"9"|) . Canonical momenta are defined by 
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The phase space (Q, P; e a , ir ea ) is equipped with canonical (nonde- 
generated) Poisson bracket. Nonvanishing brackets are {Q,P}pb — 
1, {e a , n ea }p B = 1. 

According to Eq. (1321) . in the theory there are the primary con- 
straints 7r ea = 0. Due to the condition (I30I) . Eqs. (I3T]) can be re- 
solved with respect to Q, Q — f(Q, P, e a ). Using these expressions, 
we construct complete Hamiltonian according to the standard rule 

H(Q, P, e a , 7r ea , A ea ) = 

PQ + 71 ea&a 

- L(Q,Q,e a )\\ . + A ea 7r ea , (33) 

where A ea (r) are the Lagrangian multipliers for the primary con- 
straints fl3"2"|) . Given complete Hamiltonian and the Poisson brackets, 
the temporal evolution of any quantity A is given by the equation 
A = {A,H}pb- In particular, the basic variables obey the Hamil- 
tonian equations 

Q = {Q,H} PB , P = {P,H} PB , e a = X ea , n ea = 0. (34) 

The condition of preservation in time for the primary constraints 
generally implies the secondary constraints denoted by T a 

dH 

^ea = {TTea, H}pB = = T Q = 0. (35) 

oe a 

Dynamics of the system can be equivalently obtained starting from 
the Hamiltonian action functional 



Sh — I dr 



PQ + - H(Q, P, e a , 7r ea , A ea ) . (36) 



Performing variation of the action with respect to all the variables 
(Q, P, e a , n ea , A ea ), we obtain the complete set of equations governing 
the dynamics, (EE) . AMI). (j55p . 

The procedure can be inverted. Given Hamiltonian action (136|) . 
we can restore the corresponding Lagrangian formulation. In the 
process, we can omit n ea because e a enter into the formulation 
without derivatives. Assuming that the Hamiltonian (|33|) is at 
most quadratic function of the momenta P, we write it in the form 
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H = ^PG(Q,e)P + g(Q,e). We solve the Hamiltonian equations 
Q = GP with respect to P. The solution reads P = GQ, where G 
is the inverse matrix for G. We substitute these P back into Eq. 
fl36|) . obtaining the Lagrangian action 



S= dr 



\QGQ-g(Q,e) 



(37) 



Applying the hamiltonization procedure to the action, we expect to 
arrive back at the Hamiltonian formulation ( |36|) . We do this for our 
model in Sect. 5. 

We are interested in to construct the variational problem which 
implies the constraints (I14p . ( 1T5|) and ( 128|) . The analysis made above 
suggests to take the following Hamiltonian actiorH 



S H = J dr + tt a uj a + n ei ei - H) 



(38) 



with the Hamiltonian constructed on the bas^l of constraints T 2 , T 3 
and T 4 

|(2VJ 5 " + mch) + |[(tt a ) 2 + as] + j[(co A ) 2 + a 4 ] + X £l n er (39) 

Variation of Sh with respect to momenta gives us the Hamiltonian 
equations for position variables, 

= => = ^ = e 2 (u;V - vrV), (40) 
5S 



H = o cj" = + e 2 w V, ( 41 ) 



<5S 



# _ n ^ ,-,5 



u; 5 = e 3 vr 5 + e 2 p^, (42) 



8 The Hamiltonian (139 I t is slightly different from those of the work |26| . It will lead to more 
simple Lagrangian. Both theories have the same physical sector. 

9 Peculiar property of the model is that the equation T5 = appears as the third-stage 
constraint, from the condition of preservation in time of the constraint T4. So, it is not 
necessary to include T5 into the Hamiltonian. 
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(43) 



as well as the variation of Sh with respect to the position variables 
gives the equations for the momenta, 

|| = =► p, = 0, (44) 



— ^ = => = e 2 7rV - e^, (45) 
dco fl 

vt 5 = e 2 p^ - e 4 u 5 . (46) 



Besides we obtain the primary constraints n ei = 0, which appear 
from variation with respect to the Lagrange multipliers A e; . At last, 
the variation of Sh with respect to q gives a part of the desired 
constraints of the theory, 



SS H 

Se 2 



+ mcfi = 0, (47) 



^ = =► (n A f + a 3 = 0, (48) 
be 3 

^ = =► (u/) 2 + a 4 = 0. (49) 
oe 4 



Preservation in time of the constraint (w A ) 2 + a 4 = gives the 
constraint T 5 = cj^tt" 4 = 0, which preservation, in turn, leads to the 
e 4 — — e 3 = 0. To see this, we use equations of motion (jU]) . ( 1421) . 
|5]) and flU, 



\{u A f + a 4 ]' = 2Cj a u a = 2e 3 cu A Tc A = co a tt a = 0; (50) 

(w A 7t a )- = e 3 (vr A ) 2 - e 4 (w A ) 2 = =► e 4 - — e 3 = 0. (51) 

a 4 

Time derivative of the constraint e 4 — — e 3 = determines the La- 
grangian multiplier A e4 , A e4 = fjA 3 . Preservation in time of the 
constraints T 2 = and T 3 = gives no new equations. 
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In the result, we have two pairs of second-class constraints, — 
0, 7r e4 = and T 4 = 0, T 5 = 0, while 7r e2 = 0, n e3 = 0, 

T 2 = and T 3 = T 3 + ^-T^ = represent first-class constraints. 
The first two of them are primary, and according to the general 
theory [201 12D 1221 [27] it indicates on invariance of Sh with respect 
to two-parameter group of local transformations. 

One of them is the well-known reparametrization transformation, 
its infinitesimal form is 



SY = aY, Se, 



[aei 



Sir, 



0, sx. 



(52) 



where Y = (x^p^, u A , ttb) and a = a(r) is an arbitrary function. 
Variation of Sh with respect to (|52|) is equal to the total derivative 



5S 



H 



drlaip^df + tc a uj a - -eft)]. 



The other symmetry is given by 10 l 

5x» = 0, % 



(53) 



(54) 



5u A = £lo j 



5ti j 



e 3 



(55) 



5e 2 = 0, 5e 3 = 2£e 3 , <5e 4 



e 3 



(56) 



5X ei = (5ei)\ 5n e 



0. 



(57) 



We have denoted £ = e(a 3 e3 — 0464) + 25(0363 — 0464) , where e = e(r) 
is an arbitrary function. Variation of S# is equal to a total derivative 
as well 



5S 
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dr 



— T 4 - e(a 3 e 3 
/e 3 



a 4 e4j 



(5? 



Comment. In the Berezin-Marinov model [8] the Dirac equation is 
implied by the supersymmetric gauge transformations. In our model 

10 It is Hamiltonian counterpart of the Lagrangian symmetry which will be discussed below, 
see Eqs. ll82l-|[84l). 
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the Dirac equation is associated with the £ -symmetry. So, it repre- 
sents the bosonic analogue of BM supersymmetric transformations. 

According to general theory [20j [2TJ [22], the local symmetries 
indicate on two-parametric degeneracy of solutions to equations of 
motion. Indeed, we note that A e2 (r) and A es (r) cannot be deter- 
mined neither with the system of constraints nor with the dynam- 
ical equations. As a consequence (see Eq. (H3"j) ). the variables e 2 
and e3 turn out to be the arbitrary functions as well. Since they en- 
ter into the equations for x^,u A , and it a, general solution for these 
variables contains, besides the arbitrary integration constants, the 
arbitrary functions e a . The variables with ambiguous evolution have 
no physical meaning (21]. Hence our next task is to find candidates 
for observables, which are variables with unambiguous dynamics. 
Equivalently, we can look for the gauge-invariant variables. 

The only unambiguous among the initial variables is p M , see Eq. 
(pf4l) . x M has one-parameter ambiguity due to e 2 , while co A and 
7i B have two-parameter ambiguity due to e 2 and e^. Inspection of 
equations of motion for u A and n B allows us to construct more 
quantities with one-parameter ambiguity, they turn out to be the 
angular-momentum components 



We also point out that J is invariant under the e -transformation 
055|) . The constraints (|14)) - (|15p determine the square of the angular- 
momentum tensor, see Eq. (T2"3l . Besides, they guarantee that J 5fl 
is the time-like vector 



for positive values of 03, 04. 

To proceed further, it is instructive to compare equations for 
J^ v and J 5fl with those for x M of spinless relativistic particle 
(see, for example, [2H1EH])- If we use the parametric representation 
x M (r) = (ct(r), x 1 (t)) for the trajectory x l (t), the spinless particle 
can be described by the action 



= e 2 (p fl J 5u ~p u J^) 

= e 2Pv r». 



(59) 
(60) 



(J 5M ) 2 = _ 4(a3(w 5 ) 2 + a4(7r 5 ) 2 )<0) 



(61) 




(62) 



It implies the Hamiltonian equations 



x 1 * = ep M , 



(63) 
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as well as the constraint 



p 2 + m 2 c 2 = 0. (64) 

The auxiliary variable e enters into general solution for x M (r) as an 
arbitrary function. The ambiguity reflects the freedom in the choice 
of parametrization for the particle trajectory 

x »(t)^x>»{t>) = x»{t), then l (65) 

e(r) ^e'(r') = (l + a)e(T). oe - { ae). 



The action f[62|) turns out to be invariant under the reparametriza- 
tions. 

By construction, the expression for the physical trajectory x l (t) 
is obtained resolving the equation x° = x°(t) with respect of r, 
t = t(x°), then x l {t) = x*(r(x )). The last equality implies 

d 4- = 4 = 4 (66) 

dt x° p° ( J 

As it should be, the physical coordinate x l (t) has unambiguous evo- 
lution. 

To see physical meaning of the constraint ( 1S~4"1) . we take square of 
Eq. ( 1661) and use (1S~4"1) to estimate the particle's speed 



i jf\\^-^f—^( i f] <<? . ( 67 ) 
dt ) (p 1 ) 2 + m 2 c 2 \dt ) v ' 

Hence the constraint ( 164")) guarantees that the particle's speed can 
not exceed the speed of light. 

The same result can be reproduced in the Lagrangian formula- 
tion. Indeed, variation of the action 062 p with respect to e implies 
that x M is the time-like vector, (x M ) 2 = —e 2 m 2 c 2 < 0. This also 
allows us to estimate the particle's speed 

Let us return to the spinning particle. We note both x^ and J AB 
are invariant under £ -transformations (|54p . ( I55p . So the ambiguity 
presented in equations of motion (IIP]). fl5T?|) and ([6"Uj) is due to the 



1G 



reparametrization symmetry (1521) . In accordance with this observa- 
tion, we can assume that the functions x M (r), J AB {r) represent the 
physical variables x l (t) and J AB (t) in the parametric form. Then 
equations of motion for the physical variables read 

dx % J 5i 

It = C J^' (69) 

dt x° J 50 v ; 

— ; — — c^tt — 2c — — — . (71) 
dt x° J 50 v ; 

As it should be, they are unambiguous. General solution to these 

equations will be obtained in Sect. 6. 

Although there is no the mass-shell constraint p 2 + m 2 c 2 = in 

our model, our particle's speed cannot exceed the speed of light. To 

see this, we take square of Eq. (|69i) and use the fact that J 5m is the 

time-like vector, see Eq. (1ST]) , to estimate the particle's speed 

(dx*\ 2 _ t? {J 5i y _ c2 (J 5i ) 2 



\dt J (J 50 ) 2 (J 5i ) 2 + 4(a 3 (u; 5 ) 2 + a4(7r 5 ) 2 ) 

Note that the spinning particle has causal dynamics for both positive 
and negative values of p 2 . 



4 Lagrangian formulation 

In this section we reconstruct the configuration-space formula- 
tion of the theory (1551) . We obtain various equivalent forms of the 
Lagrangian action and analyze the Lagrangian equations. While 
this is much less systematic procedure as compare with the Dirac 
method, at the end we arrive at the essentially the same results 
as those of the previous section. In what follows, we suppress the 
four-dimensional indexes, for example, we write x^u^ = (xu). We 
introduce also the following condensed notation: P a = 7r„, 7r 5 ), 
Qa _ ^x^^^^^ an d a,2 = inch. Then the Hamiltonian (1391) reads 

H = \p a G ab P b + \e A (u A ) 2 + % h 1 = 2, 3, 4, (73) 
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where G ab (cu A , e{) is a 9x9 non-singular matrix, det G = e^es^u 5 
which is schematically written as 

0(4x4) 

(r> = | e 2 w 5 r/£Tx4) e 3 77(4 X 4) 

e 2 W (lx4) °(lx4) 



A\2 



e 2 w 5 r7(Tx4) e 2^ ( 4xi) 

U(4xl) 
— e 3(lxl) 



(74) 



The notation 0(4 X 4) indicates that the first block of the matrix G ab 
is composed of the null 4 x 4-matrix, and so on. 

To find the Lagrangian, we write Hamiltonian equations for the 
position variables Q a , Q a = G ab Pb, and resolve them with respect to 
P a , P a = G a bQ b , where G is the inverse matrix for G. We substitute 
these P a back into the Hamiltonian action fl38|) . which gives the 
desired Lagrangian 



^g a g b -| e/ -i e4 (^) 2 . 



(75) 



Hence the problem of restoring the Lagrangian formulation is re- 
duced to obtaining the inverse matrix for G ab . It reads 



G, 



ab 



( G 



e3 



£3 



G 



V u e 3 (u) A ) 2UJ l J - UJl ' 



— . 
e 3 (oj A ) 2 



(76) 



where 



G 



e 3 



(e 2 u 



b\2 



It is invertible, with the inverse matrix being 



(e 2 uj 



5\2 



'■A 



'If 



5\2 



(77) 



(7f 



We substitute these expressions into Eq. ( 1751) and obtain the mani- 
fest form of our Lagrangian 



e 3 



2(e 2 u; 5 ) 2 



i; 2 + 



1 



e 2 ur 



(ixj) + 



2e,M) 2 



e 3 
e 2 ur 



iw - 



e 2 e 3 e 4 A 2 
ya 2 -ya 3 -y[(w ) + 



(79) 
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It is manifestly Poincare- invariant. The variables u 5 , e/, are scalars 
under the Poincare transformations. The remaining variables trans- 
form according to the rule 

x'v = A V + a", = A^uA (80) 

Local symmetries of the theory form the two-parameter group 
composed by the reparametrizations 

5x^ = ax 1 *, 5u A = aui A , 6e t = (ae{)', (81) 

as well as by the local transformations with the parameter e(r) 
(below we have denoted £ = e(a 3 e 3 — a^e^) + 2e{a^ez — 0464)) 

Sx" = 0, 5u A = & A , (82) 
6e 2 = 0, 6e 3 = 2£e 3 , (83) 



5e 4 = - I f j - 2e 4 £. (84) 

To confirm the symmetries, it is convenient to rewrite L in an equiv- 
alent form by rearranging its terms as follows: 



L = \G» v Dx»Dx» + ^-(u A ) 2 - ^{u A f - 1 = 2, 3, 4. (85) 

2 2e% 2 2 



G^y is the upper-left block of the matrix G a b. Besides, in Eq. (|85|) 
it has been denoted 

Dx» = - J 5 ", (86) 

I 5 ^ = — (uA^-u; 5 ^). (87) 
e 3 

Note that J 5At is £ -invariant quantity. 
The variations (I8T]) imply 

= (aG^)' — 2aG^ u , 5Dx^ = (aD/)', 
5(e 4 [(^) 2 + a 4 ]) = («e 4 [(^) 2 + a 4 ]J. 
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Using these equalities it is easy to verify that SL = (aL)'. 

Let us consider the e -transformation. We note that besides x M 
and e 2 the following Lorentz-covariant quantities 



G, v , ^, Dx», — , -HufV-uW), (89) 

u b e 3 (urr 



turn out to be invariants of the transformation fl82|) . fl83l) with an 
arbitrary £. As a consequence, the first term of the action ( 1851) 
is invariant under e -transformation. Variation of other terms of L 
under ([82]), (ESI) reads 



= I[i((^) 2 )'- 2£a 3 e 3 - 6e 4 [(oo A ) 2 + a 4 ] - 2£e 4 (ur 4 ) 2 ]. (90) 
z e 3 

Taking £e 4 = — 2e 4 £ — Xi where \ is a function to be determined, 
we cancel out the last term in fl9"0l . We are left with 

6L = l[i((u/) 2 )-+ X [{u A f + a,]] - £(a 3 e 3 - a 4 e 4 ). (91) 
z e 3 

If we further take \ = ' ^ e ^ rs ^ ^ wo terms form a total deriva- 
tive 

5L= [-^[(w A ) 2 + a 4 ]]--e(a 3 e 3 -a 4 e 4 ). (92) 
ze 3 

At last, we choose £ = e(a 3 e 3 — a 4 e 4 ) +25(a 3 e 3 — a 4 e 4 ), then the last 
term turn into a total derivative as well 

5L = l^[(u A ) 2 + a 4 ]]- [e(a 3 e 3 - a 4 e 4 ) 2 ] (93) 
ze 3 

Let us present some other possible forms of the action (1791) . 

The action implies the kinematic constraint (u; A ) 2 + a 4 = which 
is taken into account with help of the Lagrangian multiplier e 4 . 
According to classical mechanics [301 122] , we can solve the constraint 
and substitute the result back into the action, thus obtaining 

x H (iw) H ?(^) - 770.2 - — a-3, (94) 



2(e 2 u; 5 ) 2 V 04 / e 2 w 5V y 2 2 
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where now u 5 = ±\fa& + u 2 . Now, if we omit the term uj a lua in 
(1791) . it acquires the forrr0 



2(e 2 u; 5 ) 2 \ (ur 4 ) 2 ' ' ) e 2 u 5 
e 2 e 3 e 4 A 2 
— a 2 - — a 3 - —[{u ) +a 4 ] = 

^G^x" + ^(iu) - |a 2 - |a 3 - |[( W A ) 2 + a 4 ] (95) 

It is equivalent to ( 1791 since excluding the kinematic constraint we 
arrive at the expression (I94|) . As compared with (1791 . the action 
(195]) does not involve derivatives of u 5 . 

The kinetic part of the actions presented above contains the cross- 
term (xu). Starting from the action ( 185|) . we can diagonalize its 
kinetic part in an appropriately chosen variables. To achieve this, 
we introduce the auxiliary variable <t(t) and use the first-order trick 
to replace 



(Dxu) 
(u A ) 2 



2 



by a\uj A Y + 2a(Dxu). (96) 



Then the first term in (1551 reads 



-G inj Dx^Dx v - 



{Dx? {DXU)2 



(u A ) 2 



2"^ 2(e 2 w 5 ) 2 

~w^) 2 ^ Dxf + MDxuj) + a2( " A)2 ] = 
'^i^ 2 ^^^ 2 -^ 5 ^- (97) 

Now the cross-terms contained in the expression Dx^ + au^ can be 
diagonalized as follows 

Dx 1 " + ou* = x» + au". (98) 

We have introduced the new variables (the variable x 5 will appear 
below) 

~U a e 2^ 5 a 
X » = X^ UT, 

e 3 
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£5 = _ (gg) 

e 3 



e 2 u 5 \ . , e 2 u 5 



a = a+(^-\-+^-. (100) 

Comment. Hamiltonian formulation of the theory with a -variable 
implies two second-class constraints associated with this variable. 
The derivative-dependent transformation (11001) represents an exam- 
ple of conversion of the second-class constraints, see [3T] for the 
details. 

In the new variables, the last term in (197)) reads 
- a V) 2 = - (au 5 - (^> 5 - 2 

_(^+(-^)Y = -(i 5 W) 2 , (101) 

where we have replaced the variable e2 by x 5 of Eq. ( 1991) . x 5 
together with x^ form a five-dimensional space with the metric 
Vab = (-,+,+,+,-) 

Using the expressions ( 197)) - (110 II) in Eq. (185)) . the latter acquires 
the form 

e3 *\a 2 - ^a 3 - ^[(u A y + aj . (102) 



2(w 5 ) 2 2 u 2 

This almost five-dimensional form of the Lagrangian has been ob- 
tained also in . If we use the first-order trick ( 195)) to exclude the 
5" -variable, the Lagrangian reads 

£ = \GABi A i B + + 0y 2 a 2 - |a 3 - |[(^) 2 + a,U03) 



where the five- dimensional "metric" is 

u A u B 



1 /uT X 



Gab — -tf 



e 3 \ x 5 



(104) 



In contrast to G^ u , the matrix Gab has null- vector w 5 , G ab^ b = 0, 
hence it is not invertible. 
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In the rest of this Section, we analyze the Euler-Lagrange equa- 
tions which implies the action (!85|) . They read 



r c 1 

= o G^Dx^x" + -mche 2 = 0, 
de 2 2 



(105) 



|^ = =S> G^Dx^x" + J 5 ") - -(w A ) 2 - a 3 e 3 = 0, (106) 
<5e 3 e 3 



55 
(5e 4 



(w A ) 2 + a 4 = 0, then w 71 ^ = 



A,- 



(107) 



55 
5x>* 



0^[G flv Dx v ]=0, 



(108) 



55 
<5c^ 



(uDx) 



CO 



A\2 



e 3 e 3 

G iW Dx v + e 4 u; M = 0, 



e 3 



(109) 



55 
&7 



= = 



— - -u^G^Dx" 

e 3 e 3 

2 



^G^Dx"Dx u + 



(uDx) 



A\2 



—u»G„ v Dx v + e 4 w £ 
e 3 



0. 



(110) 



The equations (11 051) - (11071) do not contain the second-order deriva- 
tive and thus represent the Lagrangian constraints. 
Equation (1 108[) can be immediately integrated out 



G IIiV Dx u = = const. 
Contracting it with G afl , this can be written in the form 



111) 



Dx" = G^pv, or x" = 1^ + G^pv, or x" = y J 5/x , (112) 
if we introduce the notation 



e 2 



;ii3) 



23 



Equation f II 1 2 j) can be compared with Eq. (HOI) . The quantity J 5fi 
represents the Lagrangian counterpart of angular-momentum com- 
ponents J 5/ \ see Eq. (TIP]) . 

Using ( lllip . the Lagrangian constraint f l 1 5 f) reads (px) + \mche2 - 
0. Replacing x with help of Eq. (I112p we arrive at the classical anal- 
ogy of the Dirac equation 

p^ + mch = 0. (114) 

The Lagrangian constraint fl 1 6 [) can be rewritten in two different 
forms. First, using the expression Dx = x — I 5 , we obtain 

G^&V - a 3 e 3 - G^I 5 " - -(co A ) 2 = 0. (115) 

The constraint ( 11071) implies the identity 

G fl J^I 5 " + -(io A ) 2 = 0, (116) 

e 3 

which can be verified by direct computations. Using it in f II 1 5 [) . we 
represent the Lagrangian constraint in the form 

G ta ,x»x v - a 3 e 3 = 0. (117) 

Using the manifest form of Guv, we separate x 2 as follows 



(x^) 2 



(ux) 2 - , 



<0. (118) 



a i 

Since x M is the time-like four-vector, the particle's speed can not ex- 
ceed the speed of light. Thus the Lagrangian constraints guarantee 
causal propagation of the particle. 

Second, using Eqs. ( 11051) and ( 11111) we can exclude x M from ( 11061) . 
Then it reads 

(pi 5 ) = -{u A f + \mche 2 + a 3 e 3 . (119) 
e 3 z 

Below we present it in a more transparent form, see Eq. ( I128p . 

To simplify the equations ( I109P and f 1 11 j) for the variables u A , 
we observe two more useful consequences of Eq. (11121) . 

Using ( lllip . (I112p as well as the manifest form of G^ v , we write 

G av Dx»Dx v = p^Dx* = p^ v p v = - ej ^[p 2 - hence 

G ilv Dx^Dx v = -%u; 5 )V - (up) 2 ]. (120) 

^3 
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Contracting (11 111) with we obtain the following expression for 
(uDx) 



(ujDx) 



A\2 



— M- 

e 3 



(121) 



These equalities together with the equation ( jllip allow us to exclude 
x M from equations (I109f) . (11 lOf) for w A 



w„ e 2 5 

ujy^ 

e 3 e 3 



- e 2 



?2 



(wp) 



e 3 e 3 



p M + e 4 o; M = 0, (122) 



- e 2 



— - — (up) 
e 3 e 3 



If we introduce the notation^ 



e 2 5 



e 3 e 3 



e 3 



e 2 5 

— w p M , 

e 3 



p'" + e 4 c^ 



0. 



7T 



(wp), 

e 3 e 3 



these equations acquire the form 

tV^ = e 2 7rV t 



(123) 



(124) 



(125) 



Ti- 5 = e 2 (vrp) -e 4 w 5 . (126) 
Under the condition (11071) . the quantities n A obey the relation 



n A u A = 0, 



(127) 



which can be verified by direct computation. Besides, e^ir tta turns 
out to be proportional to —(pi 5 ) + ^(lo a ) 2 + \mcfi. Taking into 
account flllOj) . we obtain 



7r Ti A + a 3 = 0. 



:i28i 



Thus the Lagrangian constraint ( 11061) is equivalent to ( 11281) . 

Summing up, in the notation (11131) and (11241) the Lagrangian 
equations f)105j) . ( llUfip . (I108p - (lll0p can be presented in an equivalent 
form as follows 



Pf, J 5m + mch = 0. 



(129) 



12 In the next section we confirm that the quantities p 11 and n A are just the canonical 
momenta of the Hamiltonian formulation. 
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h a ha + a 3 = 0. 



(130) 



x 



-J 5 ^, (131) 



^ = e 2 7rV - e 4 u", (132) 



7T 5 = e 2 (TTp) - e A u\ (133) 
As it should be, they coincide with the Hamiltonian equations (pE0|) . 

(HSD-dSOD- 



5 Canonical analysis and quantization 

Our Lagrangian action (1851) does not involve derivatives of the 
variables e a , hence it represents an example of singular Lagrangian 
theory. The Hamiltonian formulation in this case is obtained ac- 
cording to the Dirac procedure for hamiltonization of a constrained 
system. We do it here, with the aim to confirm that the Lagrangian 
(jHSJ) and the Hamiltonian ( 13 8 p variational problems are actually 
equivalent. 

For this aim, the most convenient form of the action turns out to 
be those written in Eq. fl85|) . In this form, equations for conjugate 
momenta 

P»=q^ = G,uDx^ = G^V - I 5 "), (134) 
9L _u lx e 2 5 ~ n u_^ ea K , . 

^-d^-^~7 3 uG ^ Dx ~V,~V^ P ^ (135) 

tt 5 = |^ = - - —u^G^Dx 11 = ^ - —(up), (136) 
du 5 e 3 e 3 e 3 e 3 

can be immediately resolved with respect to velocities as follows: 

w" = e 3 7r" + e 2U y, (137) 
c0 5 = e 3 7r 5 + e 2 (wp), (138) 
G^p, + I b » = e 2 (wV - wV) = ^ J 5 ^. (139) 
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The second equality in H139[) follows with use of equations (j!37p and 
(j!38p . Conjugate momenta for the variables e\ turn out to be the 
primary constraints 

dL 

n ei = — = 0; 1 = 2,3,4. (140) 

Using the equalities (I137p - (ll40p we exclude velocities from the ex- 
pression p^ + tta^ A + ^ei^i — L + Xei^ei, and obtain the complete 
Hamiltonian 

H = j(p,J 5 " + mch) + |[(tt a ) 2 + a 3 ] + j[(u A ) 2 + a 4 ] + A ei vr e (141) 

As expected, it coincides with Eq. ( 131?p . The procedure of revealing 
the higher-stage constraints have been described in Sect. 3. We 
have obtained the following chains of constraints 

7r e2 = =s> T 2 = p^ + mch = 0. (142) 



7r e3 = 0, => T 3 = (tt 4 ) 2 + a 3 = 
7r e4 = 0, T 4 = (or 4 ) 2 + a 4 = 



=> T 5 = (cu\ A ) 



e 4 --e 3 = 0, =S> A e4 = ^A e3 . (143) 
ct 4 o 4 

The constraints T 2 and T 3 = T 3 + ^-T 4 form the first-class subset. 

If we use the Dirac bracket to take into account the constraints 
(I143p . the Hamiltonian can be presented in terms of e -invariant 
variables. First, we take into account the second-class pairs e 4 — fj-e 3 , 
7r e4 and T 4 , T5. Denoting the constraints by Ki, the corresponding 
Dirac bracket of two phase-space functions / and g is 

{/, 9h = {/> 9} - {/, K m Y\K m , g} = {f, g}+ 

^{f, (u/) 2 }{uAr A , g} - ^{/, uAr A }{(^) 2 , s}+ 

{/, e 4 - — e 3 }{vr e4 , 5-} - {/, 7r e4 }{e 4 - — e 3 , 5-}, (144) 

Q 4 Q 4 

where {, } stands for the Poisson brackets. For the basic variables 
it reads 



K, vr 5 }, = V AB - j^-y^, {u A , u B }, = 0, (145) 
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iU) 



A\2 



( W V- W V) 



J 



AB 



(146) 



Second, we impose the gauge conditions e 3 = const, u 5 = const 
for the first-class constraints 7r e 3 and T 3 , and construct the Dirac 
bracket for this set of second-class functions 

{/, 9}db = {/, g}i + {/, ^ 5 }i AfrV, <?}i - {/, 7r A 7r A }i A{u; 5 , 

{/, e 3 }l{7Te3, g}l - {/, 7Te3}l{e 3 , #}l, (147) 



where 



It implies 



A\2 



'148) 



f, ,4 Bi n AB J AC UJ C G 5B A B f-iAn\ 

{U ,7T \ DB = Cr jg^- , , W } DB = U. (149) 



{7T ,7T |db 



J 



AB 



G 5AjBC nc _ ( A ^ 



We have denoted (see also Eq. (I104p ) 

,A, ,B- 



G 



AB 



,AB 



(150) 



151] 



These formulas acquire more simple form on the constraint surface 



{U A ,7T B } DB = G 



AB 



7T 



iQ5B 



{oj a ,u b } db = 0. (152) 



Ti" 



{n A ,n B } 



J 



AB 



DB 



2(Xa 



a 3 5\A B] 

a 4 7r 5 



(153) 



The Dirac bracket of the initial variables is deformed as com- 
pared with the Poisson one. In contrast, brackets of J AB keep their 
initial form. Indeed, the constraints T a are 50(2, 3) -invariants, 
{T a , J AB } = 0. On this reason, if we compute the Dirac bracket of 
the spin-tensor components, all its extra-terms vanish, { J AB , J cd }db 
= {J AB , J CD }. Thus the transition from Poisson to the Dirac 
bracket does not modify the initial spin-tensor algebra (|ll|) . 
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By construction, the Dirac bracket of a constraint with any func- 
tion vanishes identically. So, when we are dealing with the Dirac 
bracket, the constraints can be used in all the computations as strong 
equations. In particular, we can omit them in the expression ( I14ip . 
It gives the e -invariant Hamiltonian 

Hi = j(P»J 5fl + mch) + A e2 vr e2 . (154) 

The Hamiltonian equations f l4T3]) - (j4"4"|) . f )59|) . fl60|) can now be ob- 
tained according the rule A = {A, Hi} DB instead of A = {A, H}. 

As we have shown in Sect. 2, it is consistent to describe spin- 
sector of the model by the e -invariant variables J AB instead of the 
initial coordinates u A , n B . Canonical quantization of the model 
is achieved replacing the variables x M , p u and J AB by operators 
which obey the rule [ , ] = ih{ , }, where the nonvanishing classical 
brackets are 

{x»,p v } =rf v (155) 



{ J AB , J CD } = 2( V AC J BD - V AD J BC - V BC J AD + V BD J AC ). (156) 
The operators 

j 5 m = ftr^, = nr^, (157) 

obey the necessary commutators. They act on the space of Dirac 
spinors ^(a^). The only constraint which has not been yet taken 
into account is T 2 = p M J 5fl + inch = 0. Since this is the first-class 
function, it is consistent to impose the corresponding operator on a 
state vector. This gives the Dirac equation, (T^p^ + mc)^ = 0. 

The classical equations of motion flTOj) and fITT]) imply that the 
center-of-charge coordinate x 1 experiences a complicated trembling 
motion in the theory without interaction (see also the next section). 
Besides the operator x, some other versions for the position oper- 
ator in the Dirac theory have been suggested and discussed in the 
literature [121 El H5] . In our model we can construct the following 
" center-of-mass" variable 

x" = a" + r v p v . (158) 
lp A 
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It obeys the equation 

•,, _ e2mch 

ar = ep , where e = — , (159) 

2p z 

Note also that p^ represents the mechanical momentum of the x M - 
particle. The reparametrization-invariant quantity x l (i) moves along 
the straight line, ^ = %-. We propose the variable x^ as the 
Lorentz-covariant analog of the Foldy-Wouthuysen operator. 

Using components of the spin-tensor, we can construct the Pauli- 
Lubanski vector = e^ va ^p v J a ^. It has no precession in the free 
theory, = 0, and corresponds to the Bargmann-Michel-Telegdi 
spin- vector [7]. 

In the center-of-charge instantaneous rest frame, J 50 = const, J 5 * 
= 0, it reduces to = (0, S l = e^ k pjJ k). The only J ok -part of the 
angular- momentum tensor survives in this frame. 

In the center of mass frame, p^ = (p°, 0,0,0), S is proportional 
to the three-dimensional rotation generator, S* M = (0, ~p°e* J '* Jjk). 

The equation (11581) represents the phase-space transformation 
which is not the canonical one. As a consequence, the theory became 
non-commutative. Computing the Poisson brackets, we obtain 

{x^ST} = ~ (V + i^J^p a ) , {x^p v } = rT, (160) 

{#» J 5 n = \(r)^J 5a Pa - J^ p v ). (161) 
p l 

{x", J a n = — (J^V ] - V^J^P-y), (162) 

We also present the brackets with the Pauli-Lubanski vector 

{S", S u } = 2p 2 J^ - 2p^J v ^ a p a . (163) 

S»} = -\t^J p(T - ±e^p T J a a p a , (164) 
z p 

{S», J 5a } = 2e^ aau J 5 aPu , {S^, J Q/3 } = 2p I/ e^ a J® <7 . (165) 

We point out that the second term in Eq. f 1158D has the structure 
typical for non-commutative extensions of the usual mechanics, see 

PIGHJESE5]. 
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6 Solution to the classical equations of motion 



We are interested in to solve equations of motion for the gauge- 
invariant variables x l (t), p^(t) and J AB (t), they areF*l 

~dt = c j»- (166) 

IT = 2cp ^' (167) 
*JZ = 2c t^f^, (168) 

where p M = const. Besides, the spin-tensor J AB and the time-like 
vector J 5m obeys the constraints 

p^ + mch = 0, (169) 

Jab = 8a 3 a 4 - (170) 

The differential equations written above are not polynomial. To 
avoid this difficulty, we remind that they were obtained from the 
equations for coordinates presented in parametric form 

2 

= e 2Pu J^, (172) 
= e 2 (p»J 5u - p u J 5 ^), (173) 

eliminating the ambiguity due to e 2 (T). The latter equations are 
polynomial, so their analysis represents more simple task. In the 
process, we can conventionally fix e 2 (r), since according to (I166p - 
f)170p . the physical physical coordinates we are interested in do not 
depend on a particular choice of e 2 . We take e 2 = const. After inte- 
grating the equations (117ip - (11731) . we exclude r from the resulting 
expressions thus obtaining the general solution to (I166l) - (ll68p . 

Let us start from Eq. (11721) . We compute its derivative and use 
the equations (11731) and (11691) . it gives the closed equation for J 5m 

- e 2 2 p 2 f Jfl = e 2 2 mchjf, (174) 

whose solution depends on the sign of the constant p 2 . We discuss 
the two possibilities separately. 

13 Solution to equations of motion for the initial variables ui A and ir B as well as the sub- 
sequent construction on this base the physical variables x l (t), and J AB (t) are given in the 
Appendix. 



z" = ^J 5a \ (171) 
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6.1 Spinning-particle with p 2 < 0, helical motions 

In this case, the general solution to Eq. (j!74p is given by 

777 cn 

J 5m = -^-pt* + A" cos(wr) + B 11 sin(wr), (175) 
\p\ 2 

Where \p\ = a/— p 2 , oo = e2 |p| , and A^, B^ are the integration 
constants. According to Eq. (11691) . they obey the restrictions 



Pp-Af* = 0, p^ = 0. 



(176) 



We substitute Eq. (11751) into (11711) . It gives closed equation for x^ 
which can be integrated 



z"(t) =X» + e 2 



mch 
2\pf 



p^T 



2\p\ 



A^ sm( ojt) 



1 



2\p\ 



B^ cos(wr).(177) 



By construction, only the monotonic functions x°(r) are physically 
admissible, it implies 

A = B° = 0. (178) 
We also fix the initial instant to be zero, X° = 0. Hence 

2\p\ 2 



x u = ct 



e2inch n 
2\p\ 2P 



e2fnhp° 



(179) 



Similarly, if we substitute fl!75f> into the equation J 0i = e2(p°J 5t — 
p l J 50 ), it can be integrated as well 

joi = ^0, + _±_ p o A i sin ^ r ) _ L p °B l cos(u;t), (180) 
\p\ \p\ 

E°* is the integration constant. As we know, the remaining J's are 
not independent, and can be computed according the equation 

jij = (j50yl(j5ij0j _ jSjjOiy ( lgl ) 

It reads 

"1 



,, p^YPA \p\B^Art 
J 3 = h — — z h 



mch 



—pi l A j] + ^ B H E 0j] 
\p\ mchp 



sin wr 



+ 



\p\ mchp 



cos(o;r)(182) 
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We have started our computations from the equation (I174p which 
is a consequence of (I172p . To select the subset of solutions which 
obeys the initial system fll7ip -( TT73|) . we substitute (I175p . HI 80 j) and 
(I182p into the equation (11721) . It is satisfied if 

0. (183) 



The next step is to satisfy the constraints (I169p . (11701) . The first 
one have been already taken into account, it implies (I176p . The 
constraint (I170p leads to the following restriction: 

( 2 -(A 2 + B 2 ) + ( -Mrl 2 (A x Bf = 4a 3 a 4 . (184) 



\ \p\ J \mch 
Besides, J 5fJl turns out to be the time-like vector if 

— * — * ( TftcTl \ ^ 

[Acos(wr) + -Bsin(wr)] 2 < , for any r. (185) 



/ mch\ 2 B 2 < ( mc ^ X 



It implies 



a< \m)' b< \m)- (186) 

The last step is to exclude the parameter r from the expressions 
obtained. Using Eq. (I179p . we obtain the general solution to the 
equations ( TT66|) - fTT65|) 



xHt) =X l + c^rt + -^(A i sin(wt) - B l cos(wt)), (187) 
pU 2|~ 



J 0i = ^-(A l sin(wr) - B l cos(u)t)), (188) 

\v\ 

men 

J 5 Ut) = -— p i + A 1 cos(ut) + B l sm(ut), (189) 

\vY 

jvft) = 4b [, A jI + sin(wt) - ^-p ll B j] costet), (191) 

mm |p| \p\ 



where the frequency is 



mhp° ' ^ 
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and the integration constants obey the restrictions (11841) . (11851) and 

pA = 0, pB = 0. (193) 

Let us discuss dynamics of the coordinate x l (t). 
Using the equations (11931) and (I185p . we confirm once again its 
causal motion 

f ) 2 = c w + w (icos( ^ + * sin( ^ ))2 < 

2 f b| 4 {mchf _ 2 

(p ) 2 (mfrp ) 2 |p| 2 1 J 

The curve (11871) is a helix which can be considered as a superposition 
of the rectilinear motion 

x l (t)=X i + c^rt, (195) 

and the oscillatory motion 

z\t) = rpr (A* sin(u)t) - ^ cos(u)t)), (196) 

the latter is the classical-mechanical analog of Zitterbewegung. 

Both the conjugate momentum and the Dirac equation acquire 
certain interpretation in this picture. According to Eq. (I193p . the 
Zitterbewegung oscillations occur on the plane perpendicular to p. 
This is the Dirac equation that dictates the perpendicularity. 

On this plane, the trajectory is an ellipse. To see this, we take 
the coordinate system with the origin at x % and with the axis x\ 
and X2 on the plane of the vectors A and B in such a way, that the 
xi-axis has the direction of the vector B. Then B = (B 1 ,^)) and 
A = (A 1 , A 2 ). In this coordinate system the parametric equations 
of the Zitterbewegung reads 

Xi = ——(A 1 sin(a)t) — B 1 cos(cut)), 
2\p\ 

x 2 = -L^sin^t). (197) 

2 bl 

To obtain its trajectory, we ask whether the points of the curve (I197p 
can be identified with those of a second-order line, CijXiXj = 1. 
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Short computation shows that the line is given by 




° = I -i (^ 1 ) 2 +Tb 1 ) 2 ) • (19* 

Since detC = ^^^rj > 0, the line represents an ellipse. Denoting 
its semi-axis as a and b, we can write 

The last equation together with (I186p allow us to estimate the size 
of the ellipse as follows: 

9 l9 A 2 + B 2 1 fmchV , . 

As we have seen above, canonical quantization of our model in 
the coordinates x, p and J leads to the Dirac equation. Now we 
look for the coordinates which may have a reasonable interpretation 
in the classical theory. 

Let us compute the total number of physical degrees of freedom 
in the theory. Omitting the auxiliary variables and the correspond- 
ing constraints, we have 18 phase-space variables x^, p M , u A , tta 
subject to the constraints flH]) , (EES]), fl28|) . Taking into account that 
each second-class constraint rules out one variable, whereas each 
first-class constraint rules out two variables, the number of physical 
degrees of freedom is 18 — (2 + 2 x 2) = 12. Note that this is equal 
to the number of degrees of freedom of the two-body problem. 

Further, we note that the Zitterbewegung (I196p coincides with 



jOi 

juuiuiiiait; 

(I188p . Let us choose 



the evolution of the coordinate ^ of the inner spin-space, see Eq. 



jOi i 

m = - — , f = p -- (2oi) 



J i = ^ V = 4tt - nr- (202) 

2p 0' J50 p o v > 



35 



as the new coordinates in the proble The spatial coordinates 
obey the equations 

§ = <*• d i = < 203 > 

and, according to Eqs. (I187p - (ll90p . the general solution is 

& = X i + c^rt, (204) 

fit) = -;— rfA* sin(u)t) - 5 i cosfwt)). (205) 
2|p| 

As we have seen, the trajectory of J 1 is an ellipse on the plane 
perpendicular to p. Thus, behavior of the coordinates x l and J 1 is 
similar to those of the center-of-mass and the relative position of a 
two-body system subjected to a central field. 

The Dirac electron obeys the mass-shell condition (p M ) 2 = — m 2 c 2 . 
Let us consider the subset of solutions in our model with this value 
of Then for the two-particle system with a slowly moving center- 
of-mass x\ we have p° ~ \p\ = mc. Then the frequency (11921) of the 
relative position J 1 approaches to the Compton frequency, u ~ 
Besides, the size of elliptic orbit turns out to be of the order of de 
Broglie wave-length, a/o 2 + b 2 < J^ mc - 

It would be interesting to quantize the model in the coordinates 
(120 lj) . (I202p and to compare the results with those of Dirac theory 
in the Foldy-Wouthuysen representation. 

6.2 Spinning-particle with p 2 > 0, hyperbolic motions 

The general solution to (I174p is now given by 

men 

J 5 M( T ) = — '— f- + A^e^ + B>*e-" T , (206) 

\p\ 2 

where \p\ = a/p 2 ", oj = e<2\p\. The constraint (11691) implies p^A^ = 0, 
p^B^ = 0. With this J 5fJ- (r) at hands, we immediately integrate 
equations for x M and J^ u 



We stress that we have performed a transformation on the phase-space. 
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1 B^e'^, (207) 



2\p\ 

j^{t) = + pG°^ - tf^Y^- 

-^(jfB" - tfB»)e-"\ (208) 
\P\ 

where A^, X M and S 0i are the integration constants. The com- 
ponents J 4 - 7 have been found with help of (11811) . it implies the fol- 
lowing expression for S*- 7 in terms of the integration constants 

E ij = ^\P\_r A i B j _ A j B i\ / 2 09) 
men 

We have started our computations from the equation (1174ft which 
is a consequence of ( 117211 . To select the subset of solutions which 
obeys the initial system (jTTTTl - flI73|) . we substitute (I2D511 - (120511 into 
the equation (11721) . This determines S°* in terms of other integration 
constants 

S 0i = ^]-(A°B i - B°A l ), (210) 
men 

as well as implies the restriction 

A°(pB) = B°(pA). (211) 
The equations (12091) and (I210p can be unified into the expression 

2\p\ 



men 



A [ »B v] . (212) 



The next step is to satisfy the constraints (I169p . (I170p . They 
imply the following restrictions on the integration constants 

p M A" = 0, P^ = 0, then p^ u = 0. (213) 



mch\ ( \p 



2| P | J ■ AB -{^n) K- 4 x B > - <- A ° B - BM » 1 = °»<* 214 > 

In obtaining the last equation we have used the restrictions ( 12131) . 
Besides, J 5/J turns out to be the time-like vector if 

2 + A 2 e 2u)T + 2AB + £>V 2 ^ < 0, for any r. (215) 
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Since e 2uJT increases, this inequality implies, in particular 

A^A" < 0, < 0. (216) 

The solution (12071) describes a self-accelerated motion. For the 
present case, we can not exclude the parameter r analytically. For 
sufficiently big values of r, we can neglect all the terms in the ex- 
pression for x M as compared with the third term, then 

At) - 2^4V =► = *t (217) 

It gives the following asymptotic for x l {t) 

A 1 

x\t) « X* + c—t, when t -> +oo. (218) 

Since A M is the time-like vector, see Eq. (12 16 j) . we have \v\ < c. The 
hyperbolic motions are presented also in the Frenkel theory of an 
electron, see [361. 



7 Conclusion 

In this paper we have constructed the non-Grassmann mechan- 
ical model which implies the Dirac equation. The spin degrees of 
freedom live on a seven- dimensional surface emerged in the ten- 
dimensional phase space u A , ttb equipped with the Poisson brackets 
{lo a ,ti b } = rj AB (the configuration space u A turns out to be the 
anti-de Sitter space, u A UA + a4 = 0). The surface can be covered by 
the coordinates J 5fl ,J 0l ,u 5 , where the angular-momentum compo- 
nents corresponds to the fixed value of £0(2, 3) -Casimir operator. 
In the dynamical realization this is guaranteed by the first-class 
constraint T 3 = 0, see Eqs. (JHJ) and (1231 . Due to the local sym- 
metry which corresponds to the first-class constraint, the curved 
phase space has natural structure of a fiber bundle, which allow us 
to discard the gauge non-invariant coordinate u 5 . In the result, the 
gauge- invariant angular- momentum variables J 5fl , J 0t can be taken 
as coordinates of the spin surface. Canonical quantization of the spin 
surface produces both T M and T^ u -matrices, see Eqs. (j2"5|) - (l2"T|) . 

It would be interesting to describe the underlying geometry of the 
fiber bundle (as it has been mentioned above, for the non-relativistic 
spin this is the Hopf fibration). 
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Besides the geometric constraints (jbil) and ffT5l) . which deter- 
mines the spin surface, our model implies the dynamical first-class 
constraint ( 128|) . Being imposed on a state- vector, it implies the 
Dirac equation. Although there is no the mass-shell constraint 
p 2 + m 2 c 2 = in our model, our particle's speed cannot exceed 
the speed of light. This is due to the geometric constraints which 
imply the time-like character of the four- vector J 5fJ- , see Eq. ( I6ip . 
In turn, this implies the causal dynamics of the coordinate x^, see 
Eq. dHJ. 

Analyzing the general solution (jl87p ~(??) to the classical equa- 
tions of motion (j!66p - (jl70p . we have constructed the coordinates 
(I20ip . ( 1202 j) which strongly resemble the two-body problem. If we 
assume space-time interpretation of the coordinates x l and J 1 , they 
can be identified with the center-of-mass and the relative position 
of a two-body system subjected to a central field. The Dirac equa- 
tion dictates the perpendicularity of the Zitterbewegung-plane to the 
direction of the center-of-mass motion. 

Dynamics of the subset of solutions with (p^) 2 = —m 2 c 2 is in cor- 
respondence with the dynamics of mean values of the corresponding 
operators in the Dirac theory. In particular, for the two-particle 
system with a slowly moving center-of-mass x\ the frequency f 1 1 9 2 j) 
of the relative position J 1 approaches to the Compton frequency, 
Cj ~ Besides, the size of elliptic orbit turns out to be of the 

order of de Broglie wave-length, ycfl + b 2 ~ 

Our model shows the same undesirable properties as those of 
Dirac equation in the RQM interpretation. We finish with a brief 
comment on a modification which solves the problems. We recall 
that the Dirac equation (JJ]) implies the Klein-Gordon one. In con- 
trast, in classical mechanics the corresponding constraint f[2"gj) does 
not imply the mass-shell constraint p 2 + m 2 c 2 = 0. So, the model 
presented here is not yet in complete correspondence with the Dirac 
theory. The semiclassical model that produces both constraints has 
been discussed in the recent work [23J. The extra first-class con- 
straint implies that we are dealing with a theory with one more 
local symmetry, with the constraint being a generator of the sym- 
metry [2H [221 [2Z1 EZ] • This leads to a completely different picture 
of the classical dynamics. The variable x M is not inert under the 
extra symmetry. Being gauge non-invariant, x M turns out to be an 
unobservable quantity. The variable 5 M of Eq. (I158P is gauge in- 
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variant and should be identified with the position of the particle. 
Because p^ is a mechanical momentum for the particle's speed 
cannot exceed the speed of light. In the absence of interaction it 
moves along the straight line. Hence the modified model is free of 
Zitterbewegung. 
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Appendix 

In Section 6 we have solved equations of motion for the variables 
J AB (r) and x M (r). Here we present solution to equations of motion 
for the initial variables u A , n B and x M . We shall restrict ourselves 
to the case p 2 < 0. We keep the notation p 2 = — \p\ 2 . As discussed 
before, the equations for J AB (t) and x l (t) do not depend on the 
auxiliary variables 62,3. So, we are free to select them as we want. 
We take e2,3 = 1 and define a = The equations of motion for uj a 
and ir B then read 

= tt" + u V, ^ = -aw" + 7T V; (A. 1) 

w 5 = 7r 5 + poj, 7r 5 = -au 5 + p-n. (A.2) 

Let us discuss one possible way to decouple the system above. The 
contraction of the equations for w M and 7r M with p^ gives 

(pa;)" = pn — |p| 2 Co> 5 , (jpn)' = —apco — |p| 2 vr 5 . (A. 3) 

It is suggestive to look for combinations of variables that decouple 
the system formed by (1A.2I) and ( 1A.3I) . for the variables w 5 , 7r 5 , pu 
and pn. We point out that once the equations for u 5 , 7r 5 are solved, 
one can promptly write the solutions for w' 1 and since they obey 
a closed second order differential equation 

w" + aw" = ^(tt 5 + w 5 ), (A.4) 
tt" + avr^ =^(tt 5 -au 5 ). (A.5) 
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The system (1A.2I) - (1A.3|) may be solved in a sequence of steps. First 
we define 



7T° ± = 7T° ± 



pit 

IpI 



(A.6) 
(A.7) 



The second order equations for the variables above are closed for 
the pairs (o;^,7r^) and (a/, 71"+). In fact, 



7f 5 



and 



^ = -2|p|7ri-(|p|^ + aK, 
(|p| a + a)7d - 2a|p|<4, 

u 5 _ = 2\p\ttI - (|p| 2 + a)ut, 
+ a)nl - 2a\p\uj 5 _. 



(A.8) 
(A.9) 



(A.10) 
(A.ll) 



The second step for solving the system fl A. 2[) - f)A.3[) is given by the 
rotations, 



for (lA~8l)-(fA~9l and 



y+ 
y- 



y/aUJ 5 + + 7T 5 _, 
7T 5 



+ ■ 



71 



+ ■ 



(A.12) 
(A.13) 



(A.14) 
(A.15) 



for (IA.10j) - (IA.llj) . Both systems are decoupled. The equations for 
z± and y± are 



z± = -uj±z± } 

y± = -w±y±, 



(A.16) 
(A.17) 



where u± = \p\ ± ^Ja. The general solution for z± and y± is given 
by 



z + = Acos(uj + t) + £> sin(w + r), 
Z- = C cos(w_t) + D sin(o;_r), 
y + = E cos(u + t) + Fsin(a; + r), 
y_ = Gcos(cj_t) + H sin(a;_r), 



(A. 18) 
(A.19) 
(A.20) 
(A.21) 
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where capital letters are constants of integration. Inverting the ex- 
pressions (jXl2]) - (TXl5l) and substituting flAl^ - f[OT]) provides the 
solution for u±, n±. In turn, we invert (1A.6j) - (1A.7j) to obtain the 



solutions 



[(A - E) cos(w+r) + (B — F) sin(w + r 
(C - G) cos(w_r) + (D-H) sin(o;_r 



7T ' 



-[(A + E) cos(w + r) + (B + F) sin(u; + T 
-(C + G) cos(u;_r) - (D + H) sin(u;_r 



pu 



4v^ 



[(A + £7) cos(w + r) + (B + F) sin(w + r 



pn 



+ (C + G) cos(w_r) + (D + H) sin(w_r 
\p\ 

-^-[-{A - E) cos(w+r) + -(B - F) sin(w+r 



+ 



+ 



(A.22) 
(A.23) 

(A.24) 



(C - G) cos(w_r) + (D - H) sin(a;_r)] (A.25) 



In order to solve the initial system of equations, we have taken 
one more time derivative: the variables y±, z± obey second order 
differential equations. So, it is necessary to check the consistency of 
the solutions (1A.22I) - (1A.25|) substituting them in (IA.2|) . lhs and rhs 
must coincide identically, which implies the restrictions 



A + E-- 
E-A 



B- 

- B 



F, 
- F. 



(A.26) 
(A.27) 



It follows that E = B, F = -A, G = D and H = -C. We are 
now able to solve the equations ( IA.4j) . (1A.5j) for and 7r M . Their 
general solution is given by 

w "(r) = K" cos( v / ar) + V 1 sin(v^r)- 
(A + B) , (A -B) 



4|p|Va 

(c + d; 



P^ COS(U + T) + 



£T COS(W_TJ 



(C-D) 
4|p|Va 



p^ sin(a;_| 



p^ sin(u;_r) 



(A.28) 



7r^(r) = AP sin(^r) + cos(Var) 

(A + B) u . . . (A-B) u . , 
-jr sm(u + r) + ' pr cos(a; + r) 



4|p|V« 



4|p|>/o 
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; ( + D V sin(w_r) - ^ . cos(u;_t) (A.29) 



4|p|v^ 4|_p| -^/a 

Substituting the expressions above in the initial system to check its 
consistency leads to 

iV M = -v^7P\ M" = v 7 ^^, (A.30) 
Pli K»= Pli If = 0. (A.31) 

Redefining the constants of integration 

—J— -* A < ^ B, (A.32) 
—J— ~* (A^33) 

we are left with the solutions 

w 5 (r) = —j=(A sin(w + r) + B cos(w + t) + C sin(w_r) + 

+Dcos(w_r)), (A.34) 

7r 5 (r) = Acos(w + r) — fisin(w + r) — Ccos(w_t) + 

+Dsin(a;_r), (A.35) 



w "( r ) = ^cos(v^r) + Z/sin^r)- 
cos(w+t) + p M sin(cj + r) — 



C D 

-,p^ cos(uj-t) + p M sin(cj_r) (A. 36) 



7r M (r) = y/aL* 1 cos(y/ar) - yfaK^ sm(y/ar)+ 
A B 
+ T~iP >l sin(w + r) + -r—rp^ cos(w + t) — 

\p\ \p\ 

C D 
-— r p M sin(a;_r) - —r ^ cos(w_r). (A. 37) 

Let us see what restrictions the constraints of the theory impose 
over the constants of integration. T 3 = (tt a ) 2 + a 3 = implies 

2 „„„2/ f-\ , „ r^2 • 2/ 



a 3 = aL cos (\/ar) + aK sin (\/ar) 



43 



-2aKL sin^r) cos(v^r) - {A 2 + B 2 + C 2 + D 2 ) + 
+2cos(2v^r)(AC + BD) + 2sm.(2y/ar)(AD - BC). (A.38) 



The constraint T4 

a 3 /a 4 ) 



(cu 



A\2 



+ a 4 = leads to (we remind that a 



— 03 = aL 2 sin 2 (x/ar) + aii' 2 cos 2 (\/ar) + 
+2airLsin(v^r) cos(v^r) - (A 2 + B 2 + C 2 + D ) — 
-2cos(2 v ^t)(AC + BD) - 2sm{2y/ar){AD - BC). (A.39) 



Adding the expressions (1A.38j) and (1A.39j) . we arrive at 
- 2a 3 = a(K 2 + L 2 ) - 2(A 2 + B 2 + C 2 + D 2 ). 
T 5 = oj a tt a = gives 



= cos(2^/aT 
1 



KL + -(EC - AD) 



+ sin (2 a/ot 



-(L — K ) + -MC + BD) 
2 a 



which must be zero throughout all the time, then 

aKL = 2(AD-BC), 
a{K 2 -L 2 ) =4{AC + BD). 

Now let us construct J AB and x M . By definition, 

J 5/i = 2(wV-7T 5 ^). 



Thus 



t5u 2(A 2 + B 2 - C 2 - D 2 ) u 

= -± -—= V+ 

\p\V a 

+2|p|a M sin(|p|r) + 2|p|^cos(|p|r) 



where 



a' 



\P\ 



[(A + C)L** + (B-D)Ki J 



— \{B + D)IS + (C - A)K"} 
\P\ 



(A.40) 



(A.41) 



(A.42) 
(A.43) 



(A.44) 

(A.45) 
(A.46) 
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To write J 5m we have used the solutions (1A.34|) - (1A.37I) . There is one 
more constraint to be taken into account: T2 = p^J 5 ^ + mch = 0. 
We have 



A 2 + B 2 — C 2 — D 2 = H!^. (A.47) 

J^ u can be found in the same way 

r v = 2(^/-wV) 

= Ef + Ef sin(|pir) + Ef cos(|p|r), (A.48) 

where 

Ef = 2 V ^(K"L U - K V L»), (A.49) 
Ef = — [(A - C)(K»p v - K v p») -(B + D)(L»p v - L^)](A.50) 

Ipi 

Ef = —[(A + C7)(2/ i p" - iy) + (5 + D)(K»p u - K v p^)\k.hl) 
Since x M = |«/ 5m , one has 

77? PA? 

a^(r) =X /i + ^ 7 ^r-^cos(b|r) + ^sin(|p|r). (A.52) 
2\p\ 2 

We have already used flA.47j) . The physical trajectory x l = x l {t) 
can be found with the same prescription as before: x° must be a 
monotonic function of r, then we set K° = L° = 0. So, we have (for 
simplicity, we take X°) 

n mch n 2|p| 2 . . 

map" 

Hence, 

= X i + c^rt - a* cos(wt) + i sintot), (A.54) 

where to = ^^s- We point out that the Zitterbewegung takes place 
in a plane orthogonal to the vector p l once the restrictions ( 1A.31I) 
are reduced to p ■ K — p- L — 0, leading to p ■ a — p- (3 — 0. One can 
also shows that the solution (1A.44|) for J 5fJ- (r) obeys the identity 

( J 5 "(r)) 2 + A[a 3 (u\r)) 2 + a 4 (vr 5 (r)) 2 ] = 0, (A.55) 
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that is J 5fl turns out to be the time-like vector. 

In this appendix we have solved the equations of motion for the 
initial variables (x fl ,u A ,TT B ) and confirmed the results obtained in 
Section 6: the variables J AB and experience Zitterbewegung with 
the angular frequency Co = 
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